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ABSTRACT 

X-ray luminosities of accreting T Tauri stars are observed to be systematically lower than those of 
non-accretors. There is as yet no widely accepted physical explanation for this effect, though it has 
been suggested that accretion somehow suppresses, disrupts or obscures coronal X-ray activity. Here, 
we suggest that the opposite might be the case: coronal X-rays modulate the accretion flow. We re- 
examine the X-ray luminosities of T Tauri stars in the Orion Nebula Cluster and find that not only are 
accreting stars systematically fainter, but that there is a correlation between mass accretion rate and 
stellar X-ray luminosity. We use the X-ray heated accretion disk models of Ercolano et al. to show 
that protoplanetary disk photoevaporative mass loss rates are strongly dependent on stellar X-ray 
luminosity and sufficiently high to be competitive with accretion rates. X-ray disk heating appears 
to offer a viable mechanism for modulating the gas accretion flow and could be at least partially 
responsible for the observed correlation between accretion rates and X-ray luminosities of T Tauri 
stars. 

Subject headings: stars: pre-main sequence — (stars:) planetary systems: protoplanetary disks — stars: 
formation — stars: coronae — stars: winds, outflows — X-rays: stars 



1. INTRODUCTION 

A survey of the X-ray emission from T Tauri 
stars in the Taurus- A uriga star-forming region by 
iNeuhaeuser et al.l (|1995l ) based on ROSAT observations 
revealed a significant difference amounting to a factor of 
2-3 between the X-ray luminosities of classical T Tauri 
stars (CTTS) still in the phase of active accretion, and 
weak T Tauri stars (WTTS) of the same mass that ex- 
hibit no accretion signatures. This finding was con- 
firmed in subsequent ROSAT studies of T aurus-Auriga 
and other nearby star-forming regions (e.g. lAlcala et all 



1997; Wichmann et al 



119971 : IStelzer fc Neuhauserll2001t 
Flaccomio et al.ll2003al ). in lat er XMM-Newton sur veys 
of the Taur us molecular cloud dTelleschi et al.ll2007t ) and 
NGC 2264 (|Dahm et al. 1 12001 . and by Chandra X-ray 
surveys of low-mass stars in the Orion Nebula Clus- 
ter (ONC) (iFlaccomio et alj|2003bt IStassun et ailf 



Preibisch et all l2005h . and NGC 2264 (jFlaccomio et al l 



2006). This general result might be considered surpris- 



ing, since the additional energy input through accretion 
might naively be expected to result in additional X-ray 
emission through ballist ic accretion shocks rather than in 
diminished X-ravs fe.g.lLamzin et aT1ll996l ; lLamzmll"l999l : 
iCalvet fc Gullbrindll998l)~ 

Several different explanations for the X-ray defi- 
ciency of accreting T Tauri st a rs have been pro- 
posed (e.g IFlaccomio et al.l l2003"bT: IStassun et al.l 12004 
Preibisch et all 120051: Uardine et all 120061 : iGudel et al l 



20071 : iGregory et al.ll2007l ). including: disruption of the 



magnetic corona by the encroaching accretion disk; su- 
pression of convection through the influence of accreted 
gas on stellar structure; reduction in differential rotation 
due to star-disk interactions; and obscuration of X-ray 
emission by accretion funnels. All of these explanations 
involve the suppression of X-rays as a consequence of ac- 
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cretion. We suggest here that, instead of accretion mod- 
ulating X-ray emission, X-ray emission itself might be a 
driver modulating the accretion. 

Recent results of two-dimensional photoionisation and 
radiative transfer calculations of the heating of pro- 
toplanetary disks by parent stellar coronal X-rays by 
lErcolano et al.l (|2008al ) indicate that the X-rays drive a 
disk coronal wind of significant mass flux. The estimated 
mass loss rate for a 0.75 M Q star with a rather typical 
ONC X-ray luminosity of 2 x 10 30 erg s^ 1 was of order 
10~ 8 Mq yr" 1 . Moreover, the mass loss was dominated 
by the region of the disk lying at a radial distance of 
10-20 AU. In principle, such a photoevaporation zone 
could drive off a significant fraction of the migrating gas 
from the outer disk and starve the inner disk of gas to 
replenish the accretion flow onto the star. The photoe- 
vaporation rate, and the extent to which inflowing gas is 
driven off the disk, should be correlated with the stellar 
X-ray luminosity. 

Here, we present new evidence that not only are the ac- 
cretors less X-ray luminous than non-accretors, but that 
the mass accretion rates of ONC CTTS are inversely 
correlated with X-ray luminosity, as als o suggested re- 
cently for Taurus molecular cloud stars by lTelleschi et all 
(2007). We propose that this could be a signature 
of X-ray photoevaporation-starved accretion flow: stars 
that are more X-ray bright are naturally accreting at 
lower rates because the supply of gas to the inner disk 
is diminished by photoevaporative dissipation. In £j2] 
we outline extentions to the disk model calculations of 
lErcolano et alJ (|2008af ) to investigate the influence of X- 
ray luminosity on photoevaporation rate, and in fJ3] we 
present a fresh analysis of the Ca II 8662 A accretion 
diagnostic and stellar X-ray luminosities for stars of the 
ONC ob served in the Chandr a Orion Ultra-Deep survey 
(COUP; lFeTgelson et al.|[2005l . The results are discussed 
and summarised in fj4] and ^5] 
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2. X-RAY PHOTOEVAPORATION OF DISK GAS 

lErcolano et a l. ( 2008a) have recently presented an es- 
timate of the photoevaporation rate of gas in a proto- 
planetary disk driven by X-ray heating. In order to gain 
further understanding of the dependence of the mass loss 
rate on X-ray luminosity, we have performed photoevapo- 
ra tion rate calcula t ions for identical disk parameters used 
by lErcolano et all (|2008af ). but for different values of the 
stellar X-ray luminosity. 

The basis of our X-ray photoevaporation estimate 
is a 2D photoionisation and dust radiative transfer 
model calculated usin g the MOCASSI N Monte Carlo code 
(|Ercolano et all 120031 120051 l2008bf ) of a prototypical 
T Tauri disk i rradiat ed by X-rays. The model of 
lErcolano et al. ) (2008a) assumed a 0.1-10 keV luminos- 
ity for the central pre-main sequence star of Lx = 
2 x 10 30 erg s _1 ; here we performed analogous calcula- 
tions for values of X-ray luminosity ranging from 2 x 10 29 
to 4 x 10 31 er g s" 1 , b racketing the value adopted by 
lErcolano etldl (|2008aD by factors of 10. The X-ray 
source was assumed to be an optically-thin, collision- 
dominated plasma with a temperature of logT = 7.2. 
The disk was assumed to have the initial dust and gas 
density distribution o f a 2D disk model computed by 
iD'Alessio et al.l (|1998| ) for a central star of mass 0.7 M Q , 
radius 2.5 R , and effective temperature 4000 K, and for 
an accretion rate M = 10 ~ 8 M^ yr" 1 , viscosity param- 
eter a — 0.01 (see, e.g., lPringlelll981f ). and total disk 
m ass 0.027 M^. In a refi nement over the method used 
bv lErcolano et al.l (|2008al ). the vertical gas density struc- 
ture was additionally iterated from the starting solution 
to preserve hydrost atic equilibrium ( HSE) using a version 
of the formalism of [Alexander et aTl (|2004T ) , modified to 
account for significant vertical extension of the disk gas. 
A more com plete description of the model calculations is 
presented by Ercolano et al. (20091). 

Follo wing lAlexander et all (|2004h and lErcolano et all 
(|2008al) . the photoevaporation rate through X-ray heat- 
ing at each radial grid point, R, was estimated by assum- 
ing the base of the flow is located at the height where the 
local gas temperature exceeds the escape temperature of 
the gas, defined as T es = GmnM^/kR, where M* is the 
stellar mass. To a first approximation, the mass loss rate 
per unit area is pc s , where p and c s are the gas density 
and the sound speed evaluated at the base of the flow, 
respectively. 

The total disk mass loss rate as a function of stellar 
X-ray luminosity is illustrated in Figure [TJ and the mass 
loss rate integrated over an annulus of unit width as a 
function of the annular radius, i?, is illustrated for mod- 
els corresponding to differe nt stellar X-ray lumino sities 
in Figure^ As noted by lErcolano et all (|2008aD . the 
mass loss rate peaks strongly in the 10-20 AU range — a 
zone far enough from the central star that escape tem- 
peratures become sufficiently low for the decreasing X- 
ray heating with increasing radius to effect the maximum 
mass loss. The mass loss rate is found to depend strongly 
on stellar X-ray luminosity, as expected. 

The photoevaporation rate for the canonical model 
with Lx — 2 x 10 30 erg s _1 is about a factor of ten lower 
than the non-HSE, fixe d density structure calculation of 
lErcolano et alj (2008a): this is due to the screening ca- 
pacity of the vertical evaporated gas flow that acts to 



partially shield the disk from X-rays. The HSE calcula- 
tions here will in fact overestimate this screening because 
the true gas density in outflowing unbound gas will be 
lower than that for the HSE solution in which the gas 
pressure has to support the overlying column. The HSE 
photoevaporation rate will then underestimate the true 
photoevaporation rate. Nevertheless, within the bounds 
of accuracy of the calculations, we find X-ray photoevap- 
oration rates comparable with observed accretion rates 
of CTTS. We also find a strong dependence of photoe- 
vaporation rate on Lx'- for a factor of 100 change in Lx 
the disk mass loss rate changes by a factor of ~ 50. 

3. CA II AND AN L x -M RELATIONSHIP FOR THE ONC 

Strong blue and near infrared lines of Ca II are formed 
purely in absorption in nearly all late-type stellar photo- 
spheres, with some sensitivi ty to stellar effective temper- 
ature and metallic ity (e.g. ISmith fc Drake! 119871 119901 : 
IChmielews ki 2000). In magnetically active stars, the 
lines are observed to be filled in by emission from the 
chromosphere. In the case of accreting T Tauri stars, 
the photospheric line in-filling is thought to be caused 
also by emission from regions of the upper atmosphere 
heated by magnetospheric accretion. This occurs to such 
an extent that the lines can be seen in emission, and the 
emission line strength of the infrared Ca II lines has been 
found an excellent diagnostic of mass accretion rate, M, 
( e.g. iBatalha et al 1 119961: iMuzerolle et alj 119981) . 

iFlaccomio et al.l (|2003bf ) u sed the Ca II 8662 A equ iva- 
lent width measurements of Hillenbr and et all <| 19981 ) to 
divide their Chandra sample of X-ray detected T Tauri 
stars in the Orion Nebular Cluster (ONC) into "high" 
and "low" accretion rates, and then to find differences in 
the X-ray luminosity functions of these two g roups (see 
also lStassun et al.ll2004 iPreibisch et al.ll2005h . Here, we 
re- visit the iHillenbrand et alJ (|1998h Ca II 8662 A equiv- 
alent width data with a view to making a more precise 
assessment of the mass accretion rates of ONC stars for 
which both Ca II and X-ray luminosity measurements 
exist. 

The relationship between the observed Ca II 8662 A 
flux, Fan, and mass accretion rate, M , has been stud- 
ied in detail by iMohantv et all (|2005l ) for a wide range 
of stellar mass. For CTTS, they find (their Eqn. 2) 

M = 0.71 log w (F C aii) ~ 12.66. (1) 

We adopt this relationship as the basis of our M calibra- 
tion for the ONC sam ple, together w i th the compilation 
of relevant data from iGetman et ahl (|2005l ). The latter 
includes X-ray luminosi ties from the COUP surv ey, Ca II 
equivalent widths from Hillenbr and et all (|1998h . stellar 
effective te mperatures a n d bol ometric luminosities up- 
dated from iHillenbrandl (|1997f ). and stellar ma s ses in - 
ferred from the evolutionary tracks of lSiess et al.l ([2000) . 

Since Eqn. [1] involves the flux of the Ca II 8662 A 
line, rather than the equiv alent width, we have conv erted 
the latter measurements of lHillenbrand et al.l (|1998h into 
fluxes using the predicted continuu m fluxes of the model 
atmospheres of lAllard et al.l (|2000D . These models span 
an effective temperature range of 3000-5000 K, and the 
logarithmic surface gravity range 3.5-5.5. The 8662 A 
continuum intensity for each star was obtained by inter- 
polation in the model grid for the appropriate values of 
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effective temperature and surface gravity. We derived the 
latter from the masses, tempera tures and luminosities 
tabulated bv lGetman et al.l(|2005h . Out of 502 stars with 
spectral types, effective temperatures and Ca II equiva- 
lent width measurements, we obtained M estimates for 
412. We further restricted the sample range to stars with 
mass 0.5Mq < M < 2.0M Q : being intrinsically more 
faint, stars of lower mass are prone to larger uncertain- 
ties in observed X-ray fluxes and samples are inherently 
less complete; higher mass stars have essentially radia- 
tive envelopes and likely exhibit quite different X-ray and 
disk behaviour. The stars of this sample with X-ray lu- 
minosity measurements numbered 120, and of these 45 
have mass accretion rates M > 0; the remainder have 
Ca II 8662 A either in absorption or with an equivalent 
width of zero. 

Since the final sample of 45 stars covers a significant 
mass range, and pre-main sequence stellar X-ray lumi- 
nosity depends on stellar mass, we studied the Lx~M 
relationship by comparing our derived mass accretion 
rates with deviations from the mean Lx vs M rela- 
tion d erived for the COUP sample by iPreibisch et al.l 
(|2005h . AlogLjf = logLx - log L X (M). This is illus- 
trated in Figure [3l where we find evidence for an in- 
verse correlation between M and AlogLx- the Spear- 
man's p correlation probability is 99% and the best- 
fit slope is —1.6. Similar anti-correlations are seen in 
Lx/Lboi vs M (correlation probability 99.94%, slope= 
— 1.9) and in L x /Lboi vs M/M (correlation probabil- 
ity 99.96%, slopc= —1.6). These data provide strong 
evidence that stars of a given mass with higher X-ray 
luminosit y have lower mass accre tion rates: the initial 
finding of iNeuhaeuser et al.l <fl995h that WTTS are gen- 
erally more X-ray bright than CTTS appears to be the 
manifestation of a continuous relation between X-ray lu- 
minosity and accretion rate. 

4. DISCUSSION 

In EJTJ it was noted that current ideas attempting to 
explain the lower X-ray luminosities of accreting T Tauri 
stars all seek to exploit possible suppression or disruption 
of coronal activity by accretion. Since the X-ray photo- 
evaporation rates we find in $2] are of a similar order of 
magnitude to mass accretion rates, we suggest here that 
accretion can instead be modulated by coronal X-ray lu- 
minosity. We find that the bulk of the disk mass loss 
occurs at a radial distance of 20 AU, and it is this evap- 
orative flow that can "starve" the viscous flow of gas to 
the inner disk. A recent di sk photoevaporation study by 
iGorti fc Hollenbacbl (|2009f) also finds the 10-20 AU re- 
gion to dominate the photoevaporative flow, with mass 
loss rates similar to those from our study (though in their 
model the flow is primarily driven by FUV heating; we 
defer a more complete discussion of the model differences 
to future work) . Were all inward gas flow to be curtailed 
at a radial distance of 10-20 AU, the disk inward of this 
point would essentially be cleared of gas on the viscous 
time scale. Fo r our canonical disk model, taken from 
iD'Alessio et al.l (|1998|) . the inner 10 AU contains only a 
few percent of the total disk mass and the clearing time 
is of order 10 5 yr. 

The inner disk clearing time also represents the 
timescale for which a given star must have a relatively 



high X-ray luminosity for the accretion flow to become 
completely starved of gas by photoevaporation. While 
it is possible that shorter periods of high X-ray luminos- 
ity could give rise to depleted radial "stripes" that are 
manifest in a lower accretion rate on reaching the very 
inner disk, the correlation between the two properties Lx 
and M could in such a case be erased during any given 
snapshot in time. 

The inner disk clearing time is important in the con- 
text of the large range of obseryed X- ray luminosities 
of T Tauri stars. IPreibisch et al.l (|2005l ) noted the very 
large scatter of observed X-ray luminosity in relation 
to other stellar parameters, such as stellar mass. In 
that case average deviations amount to ±0.7 dex from 
the mean relation. In order for photoevaporation to be 
at least partially responsible for the observed correla- 
tion between accretion rate and X-ray luminosity, the 
scatter in the latter must be due mostly to long-term 
intrinsic differences r ather than short-term variability. 
IPreibisch et al.l (|2005f ) ruled out stellar flaring and short- 
term stochastic variability as being responsible for the 
Lx scatter, largely by comparison of COUP observations 
with shorter Chandra ACIS ONC observations obtained 
several years earlier which indicated a median deviation 
of only a factor of ~ 2. We confirm this conclusion by 
comparison of stellar Lx values from C OUP with those 
from th e Chandra HRC ONC survey of iFlaccomio et ahl 
(20033): when sources of scatter due to different analy- 
sis procedures and different instrumental bandpasses are 
accounted for, the deviation from parity between the sur- 
veys only amounts to about a factor of 2. Nevertheless, 
the Lx scatter among T Tauri stars represents a possi- 
ble test of the relative importance of X-ray photoevap- 
oration: should future studies, perhaps over longer time 
baselines, find the scatter to be caused by relatively short 
term changes (activity cycles perhaps, that have not 
yet been seen in T Tauri stars), then photoevaporation- 
starved accretion could not give rise to the observed Lx~ 
M correlation. 

Finally, while it is suggested here that X-ray photoe- 
vaporation can drive the correlation between X-ray lumi- 
nosity Lx and mass accretion rate, and consequently the 
X-ray luminosity differences between CTTS and WTTS, 
the true explanation might involve more than one mecha- 
nism, and it is possible that both coronal disruption and 
photoevaporation play a role. 



5. SUMMARY 

New calculations of T Tauri disk photoevaporation due 
to stellar X-radiation find mass loss rates highly cor- 
related with X-ray luminosity and of a similar magni- 
tude to mass accretion rates. The results lead us to 
suggest that irradiation of protoplanetary disks by stel- 
lar X-radiation can starve the accretion flow onto the 
central star through photoevaporation. Such a process 
could explain the higher average X-ray luminosities of 
non-accreting T Tauri stars compared with those of ac- 
cretors. We find from an analysis of the Ca II 8662 A 
accretion diagnostic for T Tauri stars in Orion that the 
mass accretion rate is also correlated with stellar X-ray 
luminosity, as expected in the photoevaporation-starved 
accretion scenario. 
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Fig. 1. — X-ray heating-induced photoevaporation rate, M, from a 2D protoplanetary disk model e xpressed as a function of the model 
coronal X-ray luminosity of the central star, Lx, computed using the MOCASSIN program (see text and Ercolano ct al. 2008a for details). 
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Fig. 2. — The photoevaporation rate per unit radius, dS/dt, as a function of radial distance from the central star for different stellar 
X-ray luminosities, Lx- The evaporation rate peaks in the region 20 AU and at high values of Lx can compete with the inward viscous 
flow, starving the inner disk of the replenishing gas from the outer disk that contains the bulk (~ 90% ) of the total disk gas mass. 
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Fig. 3. — The deviation of ONC T Tauri X-ray luminosity from the mean L x vs. stellar mass, M, relation of Preibisch et al. ( 20051) as 
a function of derived mass accretion rate (see text). The Spearman's p correlation probability is 99%. The correlation indicates that, on 
average, stars with higher X-ray luminosity have lower mass accretion rate, M. 



